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Abstract 

We investigate the possibility, in nuclear fragmentation, to extract information on 
nuclear density at break-up from fragment kinetic energy spectra using a simulta¬ 
neous scenario for fragment emission. It is found that a decrease of peak centroids 
for kinetic energy spectra of fragments with increasing excitation energy can be 
observed at constant low density, which is different from recently published results 
of Viola et al. [1]. 
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One of the most challenging tasks of nuclear physics in the last decades is 
the determination of the phase diagram of excited atomic nuclei. Despite the 
important theoretical and experimental work already done, the problem is far 
from being solved. From the experimental point of view the localization of 
nuclear multifragmentation data in the phase diagram requires accurate in¬ 
dependent measurements of temperature and density at the break-up stage. 
While the problem of temperature determination has been solved with accept¬ 
able accuracy up to 5-6 MeV [2,3,4,5], no satisfactory method to determine the 
spatial extension of the presumably equilibrated nuclear system at break-up 
has been proposed. Thus, experiments using light particles interferometry [6] 
indicate freeze-out densities ranging from less than po/10 to po/2.5 (where po 
is the normal nuclear density); on the other side, statistical [7,8,9,10,11] and 
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dynamical models [12,13] succeed to describe well the available experimental 
data with freeze-out densities in the interval po/9 to po/2.5. 


Ref. [1] tries to obtain information on break-up density using kinetic en¬ 
ergy spectra of intermediate mass fragments (IMF: 3 < Z < 15) measured 
in light ion induced multifragmentation of gold, namely "^He-|-^^^Au at 50 
MeV/nucleon [14], ^"^N-|-^^'^Au at 20-100 MeV/nucleon [15] and ^He-|-^®'^Au 
at 4.8 GeV bombarding energy [16]. 


The pattern of IMF kinetic energy spectra led the authors of Ref. [1] to £t the 
extracted equilibrium sources with a Maxwell-Boltzmann type distribution. 
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where K is the kinetic energy of the considered cluster {Ap,Zp) emitted by the 
source {Ag^Zs), the cluster kinetic Coulomb energy and Vc the Coulomb 
barrier between the emitted fragments and the residual nucleus. 
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Tg is the temperature of the multifragmenting source. 

Thus, the interpretation of the behavior of these spectra with the rise of excita¬ 
tion energy is made within a parameterization suitable for sequential particle 
emission. From Eq. (1) results that a temperature increase will determine 
a shift of the centroids of the Coulomb-like peaks toward higher values of 
K together with the broadening of the distribution while a decrease of the 
Coulomb barrier (by increasing the fragments’ centre relative distances ex¬ 
pected at low density) will shift the distribution in the opposite direction. 
Starting from these premises Ref. [1] presents a systematic fit over an impor¬ 
tant collection of experimental spectra corresponding to an excitation energy 
interval ranging from 0.9 to 7.9 MeV/nucleon and reaches the conclusion that 
the displacement of the maximum of dN/dK IMF distributions towards lower 
energy and observed in the range 2-5 MeV/nucleon is a sufficient evidence in 
favor of decreasing break-up density down to ~ po/3 with increasing excitation 
energy. 

Both the short time scale characterizing the decay of nuclei with excitation 
energies exceeding 3 MeV/nucleon and the pattern of fragments’ relative ve¬ 
locities indicate that multifragmentation should be treated as a simultaneous 
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Fig. 1. Left panel: Excitation energy dependence of average values of total Coulomb 
energy, thermal energy and multiplicity of fragments with Z > 3 as a function of 
excitation energy for multifragmenting nucleus at the freeze-out density, po/S 

as obtained by MMM; Right panel: Average potential Coulomb energy experienced 
by a fragment as a function of its charge for ^®^Au, p = po/5 and Eea;=3.4 , 5.7 
and 7.9 MeV/nucleon as calculated with MMM. The Coulomb barrier experienced 
by a fragment calculated using Eq. 2 is represented with lines assuming that both 
emitted fragment and residual nucleus have normal nuclear densities (d=1.8 fm) or 
densities equal to po/5 (d=3.08 fm), see text. 

process [17,18]. In this framework, do the displacements of peak centroids of 
kinetic energy spectra reveal a decrease of the nuclear break-up density? We 
shall demonstrate in this letter that such displacements are then obtained at 
constant low density. 

To do that we shall use a microcanonical multifragmentation model (MMM) 
[11] in order to study the excitation energy dependence of the average Coulomb 
energy associated to the primary fragments at freeze-out and the IMFs kinetic 
energy spectra. To keep the treatment as intuitive as possible we assimilate 
primary fragments at break-up with spherical non-overlapping spheres placed 
in a spherical container (the freeze-out volume) and calculate Coulomb inter¬ 
action using fragment-fragment interaction, 

Ccoulomb = 1.44^^, (4) 

i<j 


where Zi denotes the charge of the fragment i, rij stands for the relative 
distance between two fragments and the sum runs over all fragments of the 
given conhguration such as to avoid double-counting. 

For simplicity we assume that for all considered cases the size of the source 
(^®^Au) and its break-up density are constant and modify only the excitation 
energy. As known from the early studies of multifragmentation, the increase 
of excitation energy induces an increase of both the degree of fragmentation 
and the thermal energy of the system. A more advanced fragmentation leads 
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to a more uniform population of the available volume and, consequently, to an 
increase of the total Coulomb energy of the system. However, by increasing 
the excitation energy, the number of fragments at freeze-out increases much 
faster than the associated total Coulomb energy which accounts for most of 
the experimentally detected hnal kinetic energy. Thus, one expects a reduced 
increase of the average Coulomb potential experienced by any fragment due to 
the mean held generated by the other fragments. These effects are illustrated 
in Fig. 1. In the left panel are plotted the average values of total Coulomb, 
thermal energies and multiplicity of fragments with Z > 3 as a function of 
excitation energy for at the freeze-out density, po/5, while the right 

panel of Fig. 1 presents the average potential Coulomb energy experienced by 
a fragment as a function of its charge, 

'ycoulomb(^) = - • 1-44 ^ 

^ iiZi=z) j Cj y{^> 


where y{Z) represents the average multiplicity of fragments with charge Z. The 
obvious relation between the total Coulomb energy Vcouiomb and the average 
Coulomb energies experienced by different fragments is, 

hcoulomb E l/('^)^Coulomb(-^)- (6) 

Z 


An increase of about 1.2 MeV/nucleon is obtained for the total Coulomb 
energy when excitation energy moves from 3.4 to 7.9 MeV/nucleon (left panel) 
whereas, as the same time, a small increase of about 0.18 MeV/nucleon is 
observed for example for Z=10 (right panel). The estimation of the Coulomb 
contribution done using Eq.(2) and d=3.08 fm which corresponds to density 
Po/5 is also shown on the right panel. This value for d is obtained taking d= 

1.8 for normal density as suggested in Ref. [1]. Estimations are indeed close 
to average values obtained considering fragment-fragment interactions. 

Adding now the kinetic part of the thermal energy at freeze-out shared at 
random between particles and fragments under constraints of conservation 
laws, we can consider what is the effect of increasing excitation energy on IMF 
average kinetic energies. The mean kinetic energy distributions as a fnnction of 
charge for the same ^®’^An source and the same density po/5 at i7ej;=3.4, 5.7 and 

7.9 MeV/nncleon are plotted in Fig. 2. At hrst glance, the behavior of < A' > 
{Z) distributions with increasing sonrce excitation is snrprising in the sense 
that while both thermal and Coulomb energies increase, the fragment average 
kinetic energies decrease. This result can be nnderstood having in mind the 
strong increase of fragment multiplicity which leads to reduced kinetic energy 
per fragment. The narrowing of < K > (Z) distribntions is obviously caused 
by the narrowing of y(Z) distribntions once the excitation energy increases. 
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Fig. 2. MMM predictions on break-up average kinetic energy as a function of frag¬ 
ment charge for source at p = po/S at £^ex=3.4, 5.7 and 7.9 MeV/nucleon. 
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Fig. 3. MMM predictions corresponding to break-up (upper panel) and asymptotic 
stage (lower panel) kinetic energy spectra for different emitted intermediate mass 
fragments resulted from the multifragmentation of ^®'^Au at po/5 and different ex¬ 
citation energies. 

Clearly these results contradict the expectation of an increase with tempera¬ 
ture or excitation energy. However they concern average quantities and not the 
peak centroids. We can consider now the spectra. As one may see in the upper 
panel of Fig. 3 the modification of the IMF kinetic energy spectra is in quali¬ 
tative agreement with the experimental data cited in Ref. [1]: with increasing 
Egx the centroids of the distribution move toward smaller energies whereas 


5 














their widths strongly increase. Since is known that primary excited fragments 
undergo secondary emission, a natural question is whether or not this process 
modihes the observed results. As one may see from the lower panel of Fig. 3 
sequential evaporations slightly diminish the IMF kinetic energies for a given 
Z, without modifying the relative displacement of distributions corresponding 
to different excitation energies. 

In conclusion, using a standard simultaneous multifragmentation model we 
explained the experimentally evidenced evolution of the IMFs kinetic energy 
spectra with increasing excitation energy as a consequence of advanced sys¬ 
tem’s fragmentation, without any assumption regarding the modihcation of 
the break-up density. To make our study as complete as possible, the behavior 
of both average kinetic energy of IMFs and the IMFs kinetic energy spec¬ 
tra have been analyzed for the freeze-out density range usually addressed by 
statistical multifragmentation models, namely po/7 to po/3. The obtained re¬ 
sults are qualitatively the same as the above presented results corresponding 
to po/5. This study suggests that an alternative explanation as compared to 
the conclusions of Ref. [1] can be proposed, which is connected to a differ¬ 
ent description of multifragmentation. Answering the important question on 
what is the break-up density dependence on the excitation energy needs more 
consideration. 


References 

[1] V. E. Viola, K. Kwiatkowski, J. B. Natowitz and S. J. Yennello, Phys. Rev. 
Lett. 93, (2004), 132701. 

[2] J. Pochodzalla et al, Phys. Rev. Lett. 75, (1995), 1040. 

[3] J. B. Natowitz et al., Phys. Rev. 65, (2002), 034618. 

[4] S. Das Gupta et al., Adv. Nucl. Phys. 26, edited by J. W. Negele and E. Vogt 
(2001) p. 89. 

[5] B. Borderie, J. Phys. G: Nucl. Part. Phys. 28, (2002) R217. 

[6] S. Fritz et al., Phy. Lett. B461, (1999), 315. 

[7] J. P. Bondorf, A. S. Botvina, A. S. Iljinov, I.N. Mishustin and K. Sneppen, 
Phys. Rep. 257, (1995), 133. 

[8] D. H. E. Gross, Rep. Prog. Phys. 53, (1990), 605. 

[9] P. D&esquelles et al., Nucl. Phys. A633, (1998), 547. 

[10] L. Beaulieu et al., Phys. Rev. G 64, (2001), 0646604. 


6 



[11] Al. H. Raduta and Ad. R. Raduta, Phys. Rev. C 55, (1997), 1344; ibid. Phys. 
Rev. C 65, (2002), 054610. 

[12] Ph. Chomaz, M. Colonna, A. Guarnera and J. Randrup, Phys. Rev. Lett. 73, 
(1994), 3512; A. Guarnera, Ph. Ghomaz, M. Colonna and J. Randrup, Phys. 
Lett. B403, (1997), 191. 

[13] J. D. Frankland et al, Nucl. Phys. A689, (2001), 940. 

[14] J. Zhang et al, Phys. Rev. C 56, (1997), 1918. 

[15] J. L. Wile et al, Phys. Rev. C 45, (1992), 2300. 

[16] D. S. Bracken et al, Phys. Rev. C 69, (2004), 034612. 

[17] D. Durand, Nucl. Phys. A630, (1997), 52c. 

[18] L. Beaulieu et al, Phys. Rev. Lett. 84, (2000), 5971. . 


7 



